The ``Spin'' Structure of the Nucleon - a lattice investigation

Preprint HLRZ 93-71 November 1993

THE "SPIN" STRUCTURE OF THE NUCLEON*
- A LATTICE INVESTIGATION

R. ALTMEYER, G. SCHIERHOLZ

DESY, Notkestrafie 85, D- 22603 Hamburg, Germany

M. GOCKELER, R. HORSLEY
HLRZ, c/o Forschungszentrum Jiilich, D-524-25 Jiilich, Germany

m

; and
0\ .

E. LAERMANN

^ Fakultat fur Physik, Universitat Bielefeld, D-33501 Bielefeld, Germany

IT} ; ABSTRACT

We will discuss here an indirect lattice evaluation of the baryon axial singlet current
matrix element. This quantity may be related to the fraction of nucleon spin carried
by the quarks. The appropriate structure function has recently been measured
(EMC experiment). As in this experiment, we hnd that the quarks do not appear
to carry a large component of the nucleon spin.

On

Oh

1. Introduction and Theoretical Discussion

Hadrons appear to be far more complicated than the (rather successful) con-
stituent quark model would suggest. For example an old result is from the nN
sigma term, which seems to give a rather large strange component to the nucleon
mass. A more recent result is from the EMC experiments which suggests that
constituent quarks are responsible for very little of the nucleon spin. These are
non-perturbative effects and so is an area where lattice calculations may be of some
help.

The EMC experiment measured deep inelastic scattering (DIS) using a /i beam
on a proton target. The new element in the experiment was that both the inital \i
and proton were longitudinally polarised. A measurement of the difference in the
cross sections for parallel and anti-parallel polarised protons enabled the structure
function g\(x. s Q 2 ) to be found. Theoretically! this is of interest as using the Wilson
operator product expansion, moments of the structure function are related to certain
matrix elements. In this case, defining s^Aq = (Ps| 57^75(7 |Ps) for q = u, d or s
quarks (i.e. the expectation value of the axial singlet current) the lowest moment
is then given by

rl . , ~ox 1/4. 1 . , 1

Talk given by R. Horsley at the 2 nd IMACS Conference on Computational Physics, Cahokia, USA.

« 0.20 + -As,

(1)

(Q 2 emc ~ HGeV 2 ). Results from neutron and hyperon decays have been used to
eliminate two of the unknowns on the RHS of this equation.

Aq can be given a physical interpretation, as the quark spin operator, Sf uark ,
is local and gauge invariant [in distinction to the gluon spin operator or orbital
angular momentum operator)tl and leads to

2(Ps\Sl uark \Ps) = AS = Au + Ad + As

« 0.68 + 3As. (2)

Before the EMC experiment, the Ellis-Jaffe sum rule was applied: the strange
content of the proton was taken as zero, i.e. As = 0. However the EMC experiment
gave for the LHS of eq. ([[]) ~ 0.126 giving As « —0.2, a large negative strange
contribution to the proton and S « 0, i.e. the total quark spin is zero. This
situation is often referred to as "the spin crisis of the quark model" , although what
was actually measured was the nucleon matrix element of the axial singlet current.

2. The lattice calculation

We now turn to our lattice calculation. By working in euclidean space on a
lattice we can turn our problem into a statistical mechanics one, which can be
approached numerically via the calculation of correlation functions using Monte
Carlo techniques. We have generated configurations^ using dynamical staggered
fermions (x) on a 16 3 x 24 lattice at j3 — 5.35, m = 0.01, which corresponds to
a quark mass of about 35MeV. Staggered fermions describe 4 degenerate quark
flavours, so we do not have quite the physical situation that we wish to describe
(i.e. 2 light quarks and one heavier quark). Nevertheless this discretisation of the
fermions has some advantages - most notably good chiral properties (as m — > 0)
when the chiral symmetry is spontaneously broken and there is a Goldstone boson,
the 7r.

The general procedure on the lattice to measure matrix elements is to consider
3-point correlation functions!. For an arbitrary operator Q, it can then be shown
that (B is the nucleon operator)

C(t;r) = <£fe*(t)n(r)fl^,(0)>

= J2 A aP (a\m^ T ^, (3)

a,(3=N,K

for \T ^> t ^> t 3> 0, where T is the time box size = 24 here, /i^ = +exp (—E N )
and the parity partner (which always occurs when using staggered fermions) is
denoted by A, so that fi\ = — exp (— E\). The amplitudes A a p and energies E a
are known from 2-point correlation functions. Thus from eq. (|3|) we can extract
(Ps\Cl\Ps). Simply setting fl to be X7i7sX nas 5 however, a number of disadvantages:

the correlation function in eq. (|3|) has connected and disconnected parts, in the fit
there are cross terms present, and the operator must be renormalised. Another
approach is to consider the divergence of the current^. This can be related to the
QCD anomaly. In the chiral limit an equivalent formulation of the problem is thus

(n f = 4)

AS = Urn \ ^ P(p),s\Q(p-p')\P(p'),s'}, (4)
p'^p i{p — p') ■ s

with topological charge

Q(p) = / d 3 xF^e^ s . (5)

(We shall take p' = 0, p = (0,0,p).) On the lattice we have used the Liischeri
construction for this charge. Although technically complicated it is integer valued
and hence has no renormalisation. In the 3-point correlation function there are no
cross terms (Q has a definite parity). A disadvantage is that we need to take p — > 0.
On our lattice the smallest momentum available is the rather large p ~ 500MeV.

Nevertheless on attempting this measurement, we find a reasonable signal with
AE rs 0.18(2). This is small and tends to support the EMC result, which would
indicate a rather large sea contribution to the proton. The existence of the QCD
anomaly also proved important. Further details of our calculation are given inS.

3. Acknowledgements

This work was supported in part by the DFG. The numerical simulations were
performed on the HLRZ Cray Y-MP in Jiilich. We wish to thank both institutions
for their support. The configurations used were those generated for the MT C project.

4. References

1. J. Ashman et. al., Nucl. Phys. B328 (1989) 1.

2. Some reviews:

S. D. Bass and A. W. Thomas, preprint, ADP-92-183/T115 (Adelaide, 1992),

H. Fritzsch, Lecture at Schladming Winter School, Schladming, 1991,

A. V. Kiselev and V. A. Petrov, preprint CERN-TH.6355/91,

E. Reya, QCD - 20 years later, Conference, Aachen , Germany, 1992, eds.

P. M. Zerwas and H. A. Kastrup, (World Scientific, 1993).

3. R. L. Jaffe and A. Manohar, Nucl. Phys. B337 (1990) 509.

4. R. Altmeyer, K. D. Born, M. Gockeler, R. Horsley, E. Laermann and G.
Schierholz, Nucl. Phys. B389 (1993) 445.

5. L. Maiani et. al., Nucl. Phys. B293 (1987) 420.

6. J. E. Mandula, Nucl. Phys. B(Proc. Suppl.) 26 (1992) 356.

7. M. Liischer, Comm. Math. Phys. 85 (1982) 39.

8. R. Altmeyer, K. D. Born, M. Gockeler, R. Horsley, E. Laermann and G.
Schierholz, Nucl. Phys. B(Proc. Suppl.) 30 (1993) 483; preprint, DESY
92-186, HLRZ 92-103, December 1992.